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Abstract: We present an NNLO-compatible subtraction scheme for computing QCD jet 
cross sections of hadron-initiated processes at NLO accuracy. The scheme is constructed 
specifically with those complications in mind, that emerge when extending the subtraction 
algorithm to next-to-next-to-leading order. It is therefore possible to embed the present 
scheme in a full NNLO computation without any modifications. 
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1. Introduction 



Exploiting the full physics potential of the LHC puts strong demands on the precise the- 
oretical understanding of QCD. In particular, accurate predictions of both signal and 
background cross sections require the computation of radiative corrections at least at the 
next-to-leading order (NLO) accuracy. In some cases however, NLO accuracy is not yet 
satisfactory and one would like to be able to compute perturbative corrections beyond 
NLO. The physical situations when this happens have been discussed extensively in the 
literature [1]. However, the straightforward (albeit toilsome) application of QCD pertur- 
bation theory to compute higher order corrections runs into the following problem. The 
finite higher order corrections are sums of several pieces which are separately infrared (IR) 
divergent in d = 4 spacetime dimensions. Therefore the evaluation of phase space integrals 
beyond LO is not straightforward because it involves IR singularities that have to be con- 
sistently treated before any numerical computation may be performed. At NLO accuracy 
the IR divergences present in the intermediate stages of calculation may be handled using 
a subtraction algorithm exploiting the fact that the structure of kinematical singularities 
of QCD matrix elements is universal. Thus it is possible to construct process-independent 
counterterms which regularize the one-loop (or virtual) corrections and real-emission phase 
space integrals simultaneously [2-5]. 

In recent years a lot of effort has been devoted to the extension of the subtraction 
method to the computation of radiative corrections at the NNLO accuracy [6-19]. How- 
ever, extending existing NLO subtraction schemes to NNLO accuracy faces the following 
problem. Consider the NNLO correction to a generic production cross section (for the sake 
of simplicity, below we will write formulae appropriate to e^e~ — > m jet production). It is 
the sum of three pieces, the doubly-real emission, the real-virtual and the doubly-virtual 
terms f r f 

^NNLO^ / da^%Jra+2+ / dcT^^^ J^+i + / da^V . (1.1) 
Jm+2 Jm+l Jm 

The above expression is formal, because da^^2 is IR singular in the one- and two-parton un- 
resolved regions of phase space, da^^^ has explicit e-poles and is singular in the one-parton 
unresolved regions of phase space, while da^ again has poles in e. Let us concentrate on 
the singly-unresolved singularities of da^^2 ™d the e-poles of duj^^^i these are just the 
divergences which are present in an NLO computation of m + 1 jet production. Thus it 
is a natural first step to try and apply a standard NLO subtraction between the first two 
terms in eq. (|1.1|) above to obtain 



^NNLO 



/ da^2 Jm+2 - d(7^2^ Jm+l + / 

Jm+2 L J Jm+l 



lm+2 

+ / daT^J 



dal\, + da,^^'^^ 



m+1 I / 



Jm+l 

(1.2) 



m 1 



where da^^2^ is just the NLO approximate cross section for m + 1 jet production and 
/i ^'^^m+2^ is ii'S integral over the phase space of one unresolved parton. 

The m + 1-parton integral in eq. ( |1.2[) above is now finite in e by construction, but this 
is of course not the end of the story yet, as it is still singular in the one-parton unresolved 
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regions of phase space. These singularities would be screened by the jet function in an NLO 
computation, but this is no longer the case at NNLO. Thus, we need to further subtract 
suitable approximate cross sections which regularize the singly-unresolved limits of dcrj^^^ 
and /-^ do"^^'^^ respectively. The construction of an approximate cross section to the real- 
virtual piece is simple enough starting from the universal coUinear and soft factorization 
formulae for 1-loop squared matrix elements [17]. But building an approximate cross 
section for da^^2^ may run into the problem that this piece need not obey universal 
factorization in the collinear limit. Indeed du^^'^^ is given by an expression of the form 

^daS5'^=da^+iC3l(e), (1.3) 

where d(T^_|_^ is the Born cross section for the production of m -|- 1 final-state partons and 
I{e) is an operator in colour space with universal pole part 

+ 0(eO), (1.4) 

the sums running over all (initial- and final-state) external QCD partons. However due to 
coherent soft gluon radiation from unresolved partons, only the sum 

(M^^UTk^M^Zi) + {M^Zi\TkTr\M^Zi) (1-5) 
factorizes in the Pi||pr collinear limit [17] {Tij. = Tj + Tr) 

C^r({M^^UTkTi\M^Zl) + {M^^U^kTrlM^Zi)) « —{M^^^\T,TMMi'^^) • (1.6) 

\ / Sir 

This factorization is violated by the factors of je^ that appear in eqs. (|l.3D and ( |l.4| ) 
at This was noticed in ref. [18] too, where it was initially claimed that the terms 

that violate factorization give a vanishing contribution after integration. However it was 
later realized that this is not the pointed out by ref. [19] first. Thus we are lead 

to define a new NLO subtraction scheme that specifically deals with this issue. We note in 
passing that it is also necessary that the unintegrated singly-unresolved approximate cross 
section have universal collinear limits, and this is not guaranteed by QCD factorization 
properties. (The approximate cross section we define in this paper will be easily seen to 
obey universal collinear factorization however.) 

For processes with only colourless partons in the initial state, such a scheme was 
presented in ref. [20], and in fact this subtraction scheme can be extended to NNLO 
accuracy [16, 17] without any modifications. 

In the present paper, we construct an NLO subtraction scheme for processes with 
hadronic initial states, that may be extended to NNLO accuracy. The steps needed to set 
up this subtraction scheme are the same as in ref. [20] . Starting from the known (tree- level) 
IR factorization formulae for singly-unresolved emission, we first write all relevant formulae 
in such a way that their overlap structure is disentangled, and multiple subtraction avoided. 



lie) 
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Then we carefully define extensions of the formulae, so that they are unambiguously defined 
away from the IR limits. As it turns out, the 'minimal' extension of the limit formulae as 
in ref. [20] does not lead to a satisfactory approximate cross section in the present case, 
and we are forced to consider a 'non-minimal' extension. Nevertheless, the structure of our 
final results will be seen to parallel those of ref. [20] closely, thus the present scheme can 
be viewed as an extension of the algorithm of that paper to processes with hadrons in the 
initial state. 

The outline of the paper is as follows. In sect.^ we recall the structure of hadronic 
NLO cross sections and the implementation of the subtraction procedure in some detail. 
In sect.m we set up some notation. Sect.Q recalls the universal factorization properties 
of QCD squared matrix elements in the collinear and soft limits. We present the new 
subtraction terms in sect.|5|. In sect.|6| we discuss the implementation of our subtraction 
algorithm and perform the integration of the approximate cross section over the phase 
space of the unresolved parton. Finally, in sect. |^ we give our conclusions. 



2. The subtraction method at NLO 
2.1 QCD jet cross sections at NLO 

Let us consider the cross section of some hadron-initiated process at NLO accuracy, 

(^{PA,Pb) = Y] [ d-qafa/AiVajfJ-p) [ ^Vb fb/ siVb, fJ-p) [^^ab^ (Pa, Pb) + (^ab^° {Pa, Pb] fJ-p)] > 

'o Jo 

(2.1) 



a,b 



where Pa = rjaPA-, Pb = ilbPB^ the subscripts on the cross sections denote the flavour of the 
incoming partons while fa/ a ^^^^ fb/B the parton densities of the two incoming hadrons, 
A and B. Assuming an m-jet quantity, the leading order parton- level cross section is given 
by the integral of the fully differential Born cross section du^^ over the available m-parton 
phase space defined by the jet function Jm 

a^^{Pa,Pb)= f da^b{Pa,Pb)Jm- (2.2) 

The NLO contribution is a sum of three pieces, the real-emission and virtual corrections 
and the collinear subtraction term 

'^ab^°{Pa,Pb;fJ'F) = da^biPa,Pb)Jm+l + da^biPa, Pb) Jm + da^biPa, Pb] tJ'p) Jm ■ 

J m+1 Jm Jm 

(2.3) 

Here the notation for the integrals indicates that the real-emission contribution involves 
m + 1 final-state partons, while the other two terms have m-parton kinematics, and the 
phase spaces are restricted by the corresponding jet functions Jn that define the physical 
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quantity. The collinear subtraction term is explicitly given by the following expression 

dcrabiPa,Pb■,^J'F) = -^Se'^ dZa / dZf, da^^^Za Pa, Zb Pb) 



X ^SbdSil - Zb) 

+ 5ac'5(l - Za 



2 \ £ 



(2.4) 



, , P'''{zb)+K^%,{zb) 



Above Sir is the phase space factor due to the integral over the {d — 3)-dimensional solid 
angle, which is included in the definition of the running coupling in the MS renormalization 
scheme* 

{27rY~^ r(l-e)' ^ ' 

The functions P°'^{z) appearing in eq. ( p.4| ) are the four- dimensional Altarelli-Parisi prob- 
abilities given by 

1 + 



X 

Pa^ix) = Tr [x^ + {l-x) 

\ \ — X 



and 
P3^[x) = 2Ca 



1 — X 



+ 1 + 3; 1 - x 



3 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



where the '+ '-distribution appearing above is defined, as usual, by its action on a generic 
test function g{x) 



dxg{x)\f{x)\. 



dx {g{x) - 5((l)]/(x) . 



(2.10) 



The factorization scheme-dependent terms Kfl^ are zero in the MS scheme. The parton 
densities are also scale/scheme dependent, so that this dependence cancels in the full 
hadronic cross section of eq. (p.l|). 

2.2 The subtraction procedure 

As is well known, the three terms of the NLO correction in eq. ( p.3| ) are all separately 
divergent, only their sum is finite for infrared-safe observables. The traditional approach 
to finding the finite NLO correction is to first continue all integrals in eq. (|2.3|) to d = 4 — 2e 



dimensions, then regularize the real-emission contribution by subtracting a suitably defined 



approximate cross section, da^^^ (j)a-,Pb)-, such that 

*The MS renormalization scheme as often employed in the literature uses Se = (47r)^e~^'''^ . It is not 
difficult to check that the two definitions lead to the same expressions in a computation at NLO accuracy. 
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• dcr^'^(pa,P6) matches the point- wise singularity structure of da^f^{pa, Pb) in the one- 
parton unresolved regions of phase space in d dimensions, 

• dcr^'^(pa,P6) can be integrated over the one-parton phase space of the unresolved 
parton analytically. 

Then, to compute the NLO cross section we write 



(^ab^'^{Pa,Pb;fJ'F) = / dal^b{Pa,Pb)Jm+l - da2'^{pa,Pb)Jn 

Jm+1 



+ / dcrabiPa,Pb) + da^biPa,Pb; + / da^f,'^(pa,pb) 

J m \_ Jl 



(2.11) 



m ■ 

e=0 



Since the first integral on the right hand side of this equation is finite in four dimensions by 
construction, it follows from the Kinoshita-Lee-Nauenberg theorem that the combination 
of terms in the second integral is finite as well for any infrared-safe observable. 

The final result is that we have rewritten the NLO contribution as a sum of two 
finite integrals, both of which are integrable in four dimensions using standard numerical 
techniques. 



Lastly, let us note that in eq. ( |2.1lD the notation da^^ {Pa,Pb)Jm is symbolic in the 



sense that da^j,' {pa,Pb) is actually a sum of terms and the jet function multiplying each 



term depends on a different set of momenta. The precise definition is given in sect. 6A 
3. Notation 

We adopt the colour- and spin-state notation of ref. [5] with a different normalization. In 
this notation the tree-level amplitude for the hadron-initiated scattering process involving 
the n final-state momenta {pi, . . . ,Pn}, 

ILyiCl,...,Cn,Ca,Cl,;Sl...,S„,Sa,Sl,/ ^ ^\ ( ^ \ 

■^n,ab [Pl, ■ ■ ■ ,Pn,Pa,Pb> , (O-lj 

is an abstract vector in colour and spin space denoted by 

\-^nib(Pi^ ■ ■ ■ ^Pn;Pa,Pb)) = (|ci, • • .,Cn,Ca,Cb) |si . . ■ , Sn, Sa, Sb)) 



(3.2) 



Here the labels a and b refer to the initial-state partons. Notice that compared to the 
definition of ref. [5] in their eq. (3.11), we do not include factors of 1/ y^njja) for each 
initial-state parton a of flavour carrying nc{fa) colours. 

Colour interactions at QCD vertices are represented by associating colour charges Tj 
with the emission of a gluon from each parton j. In the colour-state notation each vector 
\M) is a colour singlet state, so colour conservation is simply 

Y,Tj\M) = Q, (3.3) 

i 

where the sum runs over all external partons (initial- and final-state) of the vector 
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Owing to our normalization as given in eq. ( |3.2D above, the two-parton colour-connected 
squared matrix element 

• • • ,Pn■,Pa,Ph)\TiTk\M^^l^,{pl, . . . ,Pn;Pa,Pb)) (3.4) 

does not include factors of l/ndfa) (compare with eq. (3.13) of ref. [5]). The colour-charge 
algebra for the product Tj • is 

Ti-Tk=Tk-Ti, if i^k; and = Q , (3.5) 

where Ci is the quadratic Casimir operator in the representation of parton i. We have 
Cf = Tk{N^ - 1)/Nc = {N^ - l)/{2Nc) in the fundamental and Ca = 2TkNc = in the 
adjoint representation, i.e. we use the customary normalization Tr = 1/2. 

To indicate taking various kinematical limits, we use the formal operator notation 
introduced in ref. [13]. 



4. Factorization in the collinear and soft limits 

We consider the tree-level squared matrix element of a generic hadron-initiated process 
with m + 1 QCD partons in the final state \M''^\^^^^{{p};pa,Pb)\'^ ■ Here {{p};Pa,Pb) = 
{pi, . . . iPm+i'^PaiPb) denotes the set of m + 1 final-state momenta and the momenta of the 
two incoming partons. 

4.1 The collinear limit 
4.1.1 Final-final collinear 

Let us consider the limit when the momenta of the two final-state partons i and r be- 
come collinear. This limit is precisely defined by letting fej^ — > in the usual Sudakov 
parametrization of the momenta 



where is the light-like {pf^ = 0) momentum pointing in the collinear direction, while n'^ 
is an auxiliary light-like vector used to specify how the collinear direction is approached. 
The transverse momentum k'^ is orthogonal to both and {pir ■ k± = n ■ k± = 0). 
In the small-/i;j_ limit (neglecting terms less singular than l//c^), we find that the squared 
matrix element behaves as follows 

C.r\M'^Zl,abi{PhPa,Pb)\' = 

= 87^as//''^(^^^^!„,(Mm;Pa,P6)|i^^^(^^,fc±;e)|^^S,fe(Mn^;Pa,P6)) . ^^'^^ 



ir 



Here Pf.f^{zi, k±;e) are the (tree- level) Altarelli-Parisi splitting kernels (that we choose to 
label by the flavours of the two daughter partons) and the m-parton matrix elements on 
the right hand side are obtained from af,({p};Pa;P6)) by removing partons i and r 

and replacing them with a single parton ir. This parton carries the quantum numbers of 
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the pair i + r in the colhnear hmit, that is its momentum is p^^ and its other quantum 
numbers (colour, flavour) are determined according to the rule: gluon + anything gives 
gluon and quark + antiquark gives gluon. 

The Altarelli-Parisi splitting kernels are matrices that act on the spin indices of the 
parton ir. Their explicit expressions are 

,2 



{s\Pqg{z,k^]e)\s) =5ssCy- 



{s\Pgq{z,k^;e)\s) =5ssCy 



1 - z 
1 + (1 - zf 



E{l-Z) 



(/i|Pgg(z,fc_L;e)|z^) = Tr 
{^l\Pgg{z,kr:£)\v) = 2Cx 



-g^"" + 4z(l - z 



-9' 



— ez 


? 


k^k 








1 - 


z\ 



1 - z 



+ 



z 



k^ k^ 

2{l-e)z{l-z)^ 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



4.1.2 Initial-final collinear 

Now consider the case when the final-state momentum pr becomes collinear with the initial 
state momentum pa- Here the collinear limit is defined as follows 

= / I — T-lnf" -1^ h'~ i 



Par "^Pa : 



(4.7) 



1 — X 2pa ■ n 

with k± 0, and the corresponding splitting process a ^ ar + r involves a transition from 
the initial-state parton a to the initial-state parton ar together with the emission of the 
final-state parton r. The factorization formula appropriate in the collinear limit of eq. ( [4.7| ) 
reads 

Car\Ml'J^,,^,i{p};Pa,Pt)\' = 

Svras/i^" - — ({p}m ;Par,Pb)\ Pfarfr {x,k±;e)\ M''^]^^^^^^{{p}m ;Par,Pb)), 

(4.8) 



where the initial-state splitting kernels Pf^^f^{x, k_i_; e) are related to the time-like Altarelli- 
Parisi splitting kernels by 



Pf^^fA^,kr,e) = -(-l)^(/")+^(/'"-)x4^„(l/x,A;x;e) 



(4.9) 



Notice that similarly to the case of final-state splitting, we choose to label the initial-final 
splitting kernel by the fiavours of the two daughter partons. In eq. (|4.9| ) above F{q) = 1 
and F{g) = 0, which takes care of a factor of (—1) if a fermionic line is crossed. Note that 
we do not need to include the customary factor of 

w(/ar) ns{far)nc{far) 



^ifa) ns{fa)nc{fa) 



(4.10) 



in eq. ( |4.9| ) that would account for the correct counting of the number of colours and 
spins of the incoming partons, because our matrix elements are defined without these 
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factors. Nevertheless, for further reference, we record here that in conventional dimensional 
regularization 



uj{q) = ns{q)n,{q) = 2N, , and uj{g) = ns{g)nc{g) = 2(1 - e){N^ - 1) . (4.11) 

The m-parton matrix element on the right hand side of eq. ( |4.8| ) is obtained from 
l-^^m-j-i ab(iPy^Po'^Pb)) removing the final-state parton r and replacing the initial-state 
parton a with the parton ar. This parton carries the momentum par- 



4.2 The soft limit 

Suppose that the momentum of the final state gluon r becomes soft, i.e. we have = 
with fixed and A ^ 0. In this limit the squared matrix element behaves as follows 

S'-|-^2+l,a6(M;Pa,P6)P = 

where 

5.,(r) = ^ (4.13) 

is the usual eikonal factor and the m-parton matrix elements on the right hand side are 
obtained from jAI^^^^ abi{p}'^Pa^Pb)) by simply removing the soft gluon. The sums over i 
and k in eq. ( [4.12|) above run over all initial and final state partons. 

4.3 Matching the limits 

It is well known that the collinear factorization formulae contain pieces that are singular in 
the soft limit, and conversely, the soft factorization formula includes collinear divergences. 
To avoid double subtraction in the regions of phase space where the collinear and soft limits 
overlap, the common soft-collinear contributions have to be identified and the overlaps 
disentangled. A general and efficient way to do this at higher perturbative orders was 
presented in ref. [21]. However at NLO the separation of limits is quite straightforward and 
here we simply follow the approach of ref. [13] and compute the iterated limits explicitly. 

First we consider the soft limits of the collinear factorization formulae. In the case of 
the final-final collinear formula, as — > we have Zi ^ 1 and 

SrCrr\M^^\,^^,{{phPa,Pb)f =87^a,^,^' — -^T^\^^^^^ (4.14) 

Sir ^ 

if r is a gluon, while SrCir\M-^^_^^ ^f^{{p};pa,Pb)\'^ = if r is a (anti)quark. The matrix 
element on the right hand side is obtained from the original m + 1-parton matrix element 
by simply dropping parton r. The — > 0, x — > 1 soft limit of the initial-final collinear 
formula reads 

SrCar\M''^l,,^,aP};Pa,Pbt = 8na,fi^'^-^Tl\M^^^^ (4.15) 

Sar X 
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if r is a gluon, and zero otherwise. Again the matrix element on the right hand side is 
obtained from the original m + 1-parton matrix element by simply dropping parton r. 

Next we compute the collinear limits of the soft factorization formula. In the case 
where the momenta of the final-state partons i and r become collinear, p^" Zip^ and 



Pr 



(1 - Zi)p'^^, we find 



C^rSr\M^^l^^^,{{p};Pa,Pb)\'' = 8^Ta,^l 



2e 



1 2zi 



Siy 1 



(4.16) 



When the momentum of the final-state parton r becomes collinear with the momentum of 
the initial-state parton a, i.e. pr — > (1 — x)pa., we obtain 

CarSr\M^^.,A{py,Pa.p,)\''=8^a,^?^— (4.17) 



The matrix elements on the right hand sides of eqs. ( 4.16| ) and ( [4.17 ) are again obtained 
from the original m + 1-parton matrix element by dropping parton r. 

It is straightforward that the candidate subtraction term (/ and F denote the set of 
initial- and final-state partons respectively) 



Ail-M!^U.,(M;Pa,m)|' = 



(0) ^r^T.„ »,_m2 



(4.18) 



counts each singly-unresolved limit precisely once and is free of double subtractions. This 
equation is an obvious generalization of eq. (4.1) of ref. [20] to the case where coloured 
particles are allowed in the initial state. Nonetheless, eq. ( [4.18| ) above cannot yet serve 
as a true subtraction counterterm because each term in the sum is only well-defined in 
a specific collinear and/or soft limit. To define a proper counterterm, we must first give 
unambiguous meaning to all expressions away from the limits as well. 



5. NNLO-compatible subtraction terms 

5.1 Collinear subtractions 
5.1.1 Final-final collinear 



We define the extension of the collinear factorization formula of eq. (4.2) as follows 

Cf/{{p};Pa,Pb) = 

= 8na,f^''^{M'~Zbm^^^-^Pa,Pb)\PM^^^^^^^ (5.1) 

Sir 

X (1 - a,,)2<^o-2(™-i)(i--)e(ao - air-) , 
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where Pfif^ are the Altarelh-Parisi kernels as given in eqs. ( ^4.3D -( |4.q ). The momentum 
fractions Zi and Z'p are 

~ — VtQ 2^ — _ ^. — y^Q ^ (5.2) 



y{ir)Q y{ir)Q 

where we define y„Q = Suq/Q'^ = 2p„ • Q/Q^ for n = i,r and y(ir)Q = ViQ + VrQ- The 
transverse momentum fc^ ■ is 

^_|_ j Ci,rPr Cr,iPi ~l" CirPir ) Ci,r ^ ) Cr,i = -2,. '■ , (5.3) 

Oiiry(ir)Q airy(ir)Q 

where yir = Sir/Q"^ = 2pi-pr/Q'^^ while the parent momentum p^,, and air are defined below 
in eqs. (|5.5| ) and ( |5.6D respectively. Above Q is the total partonic incoming momentum in 
the laboratory center-of-mass frame, i.e. = = (^a+^6, 0, 0, K -^b)- This choice 

of transverse momentum is exactly orthogonal to the parent momentum {k±^i -pir = 0) and 
in the Pi\\pr limit behaves as required {k'j_- — > —Zi{l — Zi)sir) for any Qr- Furthermore, 
the longitudinal component of k^^i proportional to Cir vanishes by gauge-invariance when 
contracted with the matrix element, so in an NLO computation we can set C,ir = 0. However 
by choosing 

Cir- = {Zt - Zi) (5.4) 



and using the Sudakov parametrization of eq. (|4.l| ), we find that A;^ ^ ^ in the Pi| |pr limit 
(i.e. there is no residual 'gauge term' proportional to in the limit). This is necessary if 
one wants to consider a generalization of the subtraction scheme to NNLO accuracy [16]. 



The matrix elements on the right hand side of eq. (5.1) are obtained from the original 



m + 1-parton matrix element \M.^m+ii{P^'^Pa.TPh)) by removing the partons i and r and 



replacing them by a single parton ir as explained below eq. (4.2). The momenta entering 
the m-parton matrix elements, {{p}m^ \Pa7Ph) = {pi^ ■ ■ ■ ^Pir, ■ ■ ■ ,Pm+i]Pa,Pb) are defined 
as follows. First of all, it is convenient to leave the momenta of the incoming partons 
unchanged. Secondly, when defining the new final-state momenta, we find it useful to use a 
'democratic' mapping which treats the momenta of all (final-state) partons different from 
i and r identically. This choice is very convenient if one wants to extend the subtraction 
scheme to NNLO accuracy. Furthermore, in the context of an NLO computation, it leads 
to a smaller number of distinct phase space points at which subtraction terms have to be 
evaluated than the original 'dipole' mapping of Ref. [5]. We use the mapping introduced 
in ref. [20], where 

Pir = T^—ipf+Pr-a^rQ), P^ = rr^—p'^, , nj^i,r,a,b, (5.5) 



with 

_ 1 r 
air - 2 



y{ir)Q- ^Jyfir)Q-'^yir] ■ (5.6) 



Finally, note the two factors on the last line of eq. ( p.lj) . These are included to render 
the integrated counterterm m-independent (see eq. ( ^.10| ) below) and to reduce the CPU 
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time necessary for the numerical implementation [22]. For further details, we refer the 
reader to appendix A of ref. [23]. 



Phase space factorization. The momentum mapping of eq. ( |5.5| ) implements exact 
momentum conservation and leads to an exact factorization of the original m + 1 parton 
phase space of total momentum Q = pa + Pb in the form 

dcPm+ii{p};Q) = d<A^({p}(^^);Q)[dpg(p„K,;Q)] . (5.7) 

The factorized phase space measure can be written in several equivalent forms, here we 
choose the following [23, 24] 

[dpf'rl{Pr,Pir;Q)] = dai.(l - ai.)2(™-l)(l-^)-l f|^d(/<2 (pi , ; )0(«-)e(l - a^,) , 

(5.8) 

where s—^ = 2pir ■ Q and p'^.^^ = (1 - air)pfr + "irQ^- 



Integral of the subtraction term. Finally we compute the integral of the subtraction 
term over the unresolved phase space [dpi^jl{pr,Pir]Q)]- Firstly, since k'^- as given by 
eq. ( |5.3D is orthogonal to pf^., the spin correlations generally present in eq. ( |5.1| ) vanish 
after azimuthal integration [5]. Thus when computing the integral of the subtraction 
term C^^{{p};pa,Pb) over the unresolved phase space, we can replace the Altarelli-Parisi 
splitting kernels Pf^f,. by their azimuthally averaged counterparts Pf^fr- Then we have 



[dpi ;m {Pr , Pir ; Q)]Cf/ {{p] ■,Pa,Pb) 

as 



g5.^^)'cf/(y~Q;e)Ti|A^;^;„,({pe);p„p,) 

where y—^ = 2pir ■ Q/Q'^ and 
(4^)2 



,(0) 



(5.9) 



Ss 



-iQ 



2\e 



da,v (1 



' 27r 



<i<p2 {Pi ,Pr;P(ir)) — Pf,frizi,^)7f2 



(5.10) 



As claimed, C[^^ {yj^ e) is independent of m and in order to lighten the notation through- 
out the paper, we do not explicitly indicate its dependence on oq and do- 

The integral appearing in eq. (5.10) has been computed in refs. [23,24] (see also [20] 
for the case of ckq = 1) and here we only recall the final results. The pole parts of the 
integral are independent of ao and do and we have 



Cf/(x;e) 



19 



^ + |-'Mx)+0(e») 



19 



7 \ — 



2 11 

Ye 



ln(x) + 0(e° 



(5.11) 
(5.12) 
(5.13) 



- 12 - 



We present the finite parts for tlie specific choice of = 3 — 3e in appendix A.l 



5.1.2 Initial-final collinear 

The extension of the collinear factorization formula of eq. ( |4.8| ) is written as 
Ciri{p}'^Pa,Pb) = Svras/"^^-^ — 

X 0(yo - Var) , 

where Pf^,.f,. are the initial-state Altarelli-Parisi kernels as given in eq. (^^). The momen- 
tum fraction Xa is given by 

Xa = l-yrQ, (5.15) 
while the transverse momentum ^ reads 

k'l,a=P'f-yrQP^a-yarQ''. (5.16) 

Again Q is the total incoming partonic momentum in the laboratory center of mass frame, 
that is = Pa + ^^^^ thus = 2pa ■ Pb = '^Pa • Q- The choice of k'^ ^ in eq. ( 5.1(^ ) 



is orthogonal to the parent momentum par (^±,a • Par = 0) to be defined in eq. ( 5.17| ) 
below and has the correct behaviour in the Pa||Pr limit: k'j_^ —>■ —(1 — x)saT- Furthermore, 
using the appropriate Sudakov-parametrization, eq. ([4.7]), we find that in the collinear limit 
k'^a ^ ^ (again there is no 'gauge term' in the limit). This is important if the present 
NLO algorithm is to be included as part of an NNLO subtraction scheme. 



Finally, the matrix elements on the right hand side of eq. (5.14) are obtained from the 
original m -\- 1-parton matrix element on the left by removing the final-state parton r and 
replacing the initial-state parton a with the parton ar as discussed below eq. ( [4.111 ). The 
m momenta {{p}ln^;Par,Pb) = iPi, ■ ■ ■ , Pm+V, Par , Pb) (Pr is missing from the set) entering 
the factorized matrix elements are defined as follows. Firstly, it is convenient to leave the 
momentum of parton b unchanged. Secondly, it is also convenient to define the momentum 
of parton ar, par to be parallel with pa- Thus we set 

P^r = XaP^, pf^ = A[Q -pr,Q - [1 - Xa)pa]uPn^ n a,b,r , (5.17) 
where A{K, K)u is 

K{K, Kt - < ^K + kf ■ ^ ^ ^ 

This is the matrix of a (proper) Lorentz transformation, that maps k^ to K^, whenever 
= k^. This mapping treats all final-state momenta except r identically. 
The 0-function is included on the last line of eq. ( 5.14[ ) to reduce CPU time in a 

numerical implementation. 
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Phase space convolution. The analogue of the phase space factorization in eq. ( ^.T] ) is 
the following phase space convolution 



d(prn+li{p};Pa+Pb) = / dx d(l)rn{{p}'^'~'^ ] XPa + PbMp'i^miPr, Pa, x)] , 

Jo 



(5.19) 



where the one-particle phase space measure can be written as 



[dp[^%r,Pa,x)] = -^(QY--dyardia&{XaMl-Xa)&{yar)e [ 1 



l-Xa 



X (1 — x) 



-2e 



Var 



f Xa \ 1 Xq^ 



(5.20) 



5{x - Xa) . 



Above we have yar = "^Pa • Pr/Q"^- 

We note that the momentum mapping in eq. ( 5.17 ) and consequently the phase space 
convolution in eq. ( ^.19 ) coincides exactly with the mapping and convolution of section 5.5 
of Ref. [5]. Our Xa is exactly their xi^ah, after the replacement i ^ r and similarly, their Vi 
is just our yar after the same replacement. 



Integral of the subtraction term. When computing the integral of Cl^{{p};pa,Pb) 
over the one parton factorized phase space [dp^i^^^iPriPa, x)], we may replace the splitting 
functions Pfarfr with their azimuthally averaged counterparts Pf^rfr since /c^ ^ is orthogonal 
to Par-, and the spin correlations generally present in eq. ( |5.14| ) vanish after integration. We 
find 



[dpS ' P'' ' ^)]'^ar {{P] ;Pa,Pb) 



IF I 



(2 \ £ 



(5.21) 



with 



C 



a{ar) 



{x]e) 



, , . I dyar / dXa 6 ( 1 
'^(Ja) Jo Jo 



I- Xn 



(1-x) 



-2e 



Va 



1- X„ 



1 



Va 



I- Xr. 



1 1 

6{X - Xa) — Pf,,fAXa;e)-^ , 
i>ar 1 nr 



(5.22) 



where we have introduced factors of io{fa) and io{far) into the definition of C^^^^-j (x; e) above 
for later convenience. As before, to lighten the notation, we do not explicitly indicate the 
dependence of C^(^^)(x;e) on y^. 

The integral which appears in eq. ( 5.22| ) above has been computed in ref. [5] for y^ = 1 
and in ref. [25] for general y^, so here we only recall the final results. The pole part of the 
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integral is independent of and reads 



ICf 



[x' + (1 



1 fl + x' 



e \ \ — X 

2 r/ 1 

e 



1 — x 



+ 0(e°) , 

-x)2]+0(eO), 

+ (5(1 -X 
1 — X 



2e 



(5.23) 
(5.24) 
(5.25) 



+ 



1 + x(l — x) 



+ (5(l-x)^ + 0(e°), (5.26) 



where the '+ '-distribution appearing above is defined in eq. (|1|). Note that for = 1, 
these functions are in fact identical to the V"'"*(x;e) functions introduced in ref. [5], up to 
a factor of (after the usual replacement of i — > r). Explicitly we have C^j-^^-j(x; e; = 
l)Tl^ = V"''"^(x;e). The finite parts of these expressions are given in appendix |A.1|. 



5.2 Soft-type subtractions 

Next we define the extension of the soft (eq. ( |4.12[ )) and soft-collinear (eqs. ( [4.16 ) and ( [4.17| )) 
limits. We call these 'soft-type' subtractions because they all use the same 'soft' momentum 
mapping, see eq. ( 5.30 ) below. We set 

i k^i 



><(-^Sa6(M^;:^;i'a,P.)|T.T,|A^2!,,(W(^);Pa,P6)) 

X (1 - y.Q)'^o-('"-i)(i^^)e(yo - yrQ)P{yrQ) , 



Cir ■S^-i{p}',Pa,Pb) = SnasfJ, 



2e 



1 2£; 



rTf\M^^U{pVr:^-^Pa,Pb 



Sir 1 



(r).. 



Cir'Sr{{p}'^Pa,Pb) = Svras/i 



X (1 - y.Q)<>-(™-i)(i-^)G(yo - yrQ)P{yrQ) 
1 2 „o , . ,(0) 



2e 



■Tl\MZbm^^;Pa,Pb 



Sar 1 Xa 

X (1 - y.Q)"o-(— i)(i--)G(yo - yrQ)P{yrQ) . 



(5.27) 



(5.28) 



(5.29) 



The eikonal factor, Sik{r), and the momentum fractions Zi and Xa which appear in the 
above equations were defined in eqs. ( [4.13 , [5.2| ) and (5.15) respectively. (Also recall that 
the sums over i and k in eq. ( 5.27[ ) run over all initial and final state partons.) The matrix 
elements on the right hand sides are obtained form the original m+l-parton matrix element 
by simply dropping the parton r. The momenta {{p}m;Pa,Pb) = {pi, ■ ■ ■ ,Pm+i;Pa,Pb) iPr 
is missing) entering the m-parton matrix elements are defined as follows: 



p':, = A[Q,{Q-Pr)/Xr]^{p''/Xr), Xr = ^/l^^ , n^r,a,b.. 



(5.30) 



where A is defined in eq. ([5.18[) . Notice that the momenta of the incoming partons are left 
unchanged. This mapping treats all final-state momenta except pr in an identical fashion. 
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In eqs. (|537| )-(|5^) above, we have included some additional factors. The first two 



of these serve the same purpose as in the collinear case, see eq. (5.1). Firstly, they render 
the integrated counterterms m-independent (see eqs. ( ^.35 , 5.39 , 5.4C1| ) below) and secondly 
they reduce the necessary CPU time in a numerical implementation. 

The final factor of P{yrQ) is included to control the exact form of the integrated soft- 
type subtraction terms. A similar approach, albeit in a different context, has been used in 
ref. [26]. In the present case this is crucial for obtaining an insertion operator that obeys 
universal collinear factorization. In order that the terms in eqs. ( |5.27| )-( |5.29| ) above reduce 
to their correct forms in the soft limit, we must have P{yrQ) ^ 1 as y^Q 0. Hence we 
set 



TV 



(5.31) 



n=l 



where we shall find below that the minimal number of terms necessary is = 2 and we will 
also determine ai and a2- The a„ will be seen to depend on the parameters yo and as 
well as e. In order to keep the notation streamlined, we do not indicate these dependences. 



Phase space factorization. The momentum mapping defined in eq. ( 5.30| ) above leads 
to an exact factorization of the original m + 1-parton phase space of total momentum 
Q = Pa + Pfe in the form 



dcPm+ii{p};Q) = d</.„({p}W;Q)[dpi;i(p.;Q)] . 
We choose to write the factorized phase space measure as in refs. [23, 24] 



(5.32) 



[dp^;l,{pr; Q)] = dyrQ (1 - y.Q)("^-^)(i-^)-i|^d02(Pr, K; Q)G(y.Q)G(l - y.g) , (5.33) 

where the time-like momentum K is massive with = (1 — yrQ)Q'^ ■ 

Integral of the subtraction term. Consider first the integration of the soft counterterm 
'5r({p};Pa,Pb) over the one-parton unresolved phase space [dvVm{Pr:Q)]- We find 



2 \ £ 



j [dPi;i(Pr;Q)]5,,(M;Pa,m) = 
.ieF keF 

+ 2Y,Sir{Y^lQ;e){M'^^l,ml:^-,Pa,Pb)\TaT.^^^^^^ (5.34) 



i&F 



+ 2Si^(e)(>f^;,,({p}^^Pa,p,)|raT,|>f^;,,({p}^^Pa,p,)) 



<(0) 
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where a and h are the labels of the two incoming partons and 



--^(Q^y J dy,Q(l-y,Q)'^o-ip(y^Q)^ J dMPr,K;Q)-Sik{r). 

(5.35) 

Here Ji,Jk = I,F denote whether i and k are initial- or final-state partons. Note that 
^ifc"'* (^fc Q' ^) symmetric in i and k, so the ordering of indices does not matter. In 
particular S^/{Y;£) = Sl^{Y;e). As promised, S^"^''(l^^ g; e) is independent of m and 
depends on i and k only through the combination 



Y- 



y'lk 



(5.36) 



ViQy'kQ 

Notice that Pa = Pa, Pb = Pb (the soft mapping of eq. (5.30) leaves the incoming momenta 
unchanged), so UaQ = Viq = 1 and thus Y^j^q = 1. Then in eq. ( ^.34 ), S^^ only depends 
on e, as indicated. As usual, we do not show the dependence on yo and to lighten the 
notation. 

Next, consider the integral of the soft-collinear subtraction terms. We have 

,2 \ £ 



(Pr ; Q)]Ci/ -Sr {{p} ;Pa,Pb) 



FF , 



27r' 



Q 



CS^^{e)T^\MZbmt^;Pa,Pbr 



(5.37) 



and 



I [dpi;i(p,;Q)]Cff5,({p};p„,p,) = (^) CS'^{s)Tl\MtU{p}t^-,Pa,Pb)\' 

(5.38) 



with 



yo 



27r 



>{pr,K;Q) 



_]_JZi 
Sir 1 



~ 5 

Zi 



(5.39) 



and 



CS^^(e) = ^(Q')'£"dy.Q(l-2/.Q)'^o-ip(y.Q)|^ J dMPr,K;Q) 



1 - Xi 



(5.40) 

As before, CS^^(e) and CS^^{e) are independent of m but do depend on yo and d^. The 
dependence on the later two parameters however is suppressed in the notation. 

The soft integral in eq. ( ^.35 ) for i and k both final-state partons, as well as the soft- 
collinear integral in eq. (|]3|) ab ove have been evaluated in refs. [23,24], with P{yrQ) = 1 
(see also [20] for the case of yo = 1). Computing the soft integral for the cases when either 
one or both of i and k are initial-state momenta and the soft-collinear integral in eq. ( 5.4C| ) 
above, with P[yrQ) as given in eq. (5.31), is straightforward with the same techniques. We 
nevertheless present the details here, as they are needed to derive the a„ that enter P{yrQ) 
in eq. ( ^.3lD . 

First, in order to write the factorized phase space measure in eq. ( 5.33| ) explicitly, 
we choose a specific Lorentz-frame. The choice of a convenient frame is different for 
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S^^(>ifc,Q;e)> SidYalQ-^^) and S^^(e) as well as for CS^^(e) and CS^^(e). In particu- 
lar, in each frame we have 

Q'^ = ^/s{l, . . .) , and = £',.(1, ..'angles'. ., sin sin sin i? cos cos-!?) , (5.41) 

where the dots stand for vanishing components, while the notation ..'angles'., denotes the 
dependence of pr on the d — 3 angular variables that can be trivially integrated. Then, for 
computing the various integrals, we set 







= Ei{l,...,l), 


~u 
Pk 


= Ek{l,. 


..,sinx-fc,cosXifc 






= f(i.....i). 


~n 
Pi 




.,sinx„-,cos Xai) 




P'a 




Pb 


= ^(1, 
2 ^ ' 


...,-1). 



(5.44) 



Now in terms of the scaled energy-like variable 

2pr • Q 2Er 



g2 ^ 

and the angular variables ?? and ip the two-particle phase space reads 

X d(cost?) d(cos(/?)(sin7?)~^'^(sin(/?)~^~^^ , 



(5.45) 



(5.46) 



for all three choices of frame in eqs. ( |5.42| )-( |5.44| ). In eq. ( ^.46| ), the limits of integration on 
Er are and 1, while the cosines of both angles run from —1 to 

To write the integrands in these variables, we observe that the precise definitions of 
and Pk as given in eq. ( 5.30| ) above imply 



^ik — (1 ^^)^ik ' — ^ir ' ^f'^ — ^kr 



Si, 



iQ 



(1 - er)sjo + s-„ , (5.47) 



and 



iQ 



(5.48) 



Then in the Lorentz-frame of eq. (|5.42|) we find 



Sik _ ^^ ^^ s-.f. _ 4:Y-.-^ Q (1 - 



^ir^kr 

and 



Q'^ (1 — cos t?)(l — cosx^j^ cos — sinx^^ sini? cos 99) ' 

(5.49) 



1 Zj 

Sir 1 



Sir ■SrQ 



1 1 



1 + 



2(1 -g.) 

er(l — cos ■!?) 



(5.50) 
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Similarly, in the frame of eq. ( |5.43| ) we obtain 

1 SaQ 



1 



1 



/-y ^ ai _j 

'^KiQ yr^r 1 

el (1 — COST?) (1 — cos Xfli cos 1? — sin sin cos 



+ 



1 



\/l - Er 



EriX + \/l — £t)Q'^ \ 1 ~ COS I? 1 — COS Xai ^'^^ ~ Xai ^ ^ 
1 



(1 + ^/r^)2Q2 



1 1 



11 2 1 1 



Sar 1 - ia Sar VrQ 1 - COS 1? 



and 



Finally, in the frame of eq. ( |5.44| ) we have simply 

Sab 4 1 1 



SarSbr Q"^ E^ 1 — COS^ 1? 



(5.51) 
(5.52) 

(5.53) 



The Y--j^ Q and Y^^q that appear in eqs. ( 5.4g| ) and ( 5.51 ) above are defined in eq. ( 5.36 ). 

Now we are ready to compute the integrals in eqs. ( p. 35 , 5.39) and ( |5.40| ). We start 
by computing S^^^(%Q;e) and CS^^(e). Using eqs. (|5l|, [slol) and (|5^ ) we obtain 



N 



By^ (-2e, d'o + 1) - ^ a„ By, (n - 2e, d'o + 1) 



(5.54) 



and 



n=l 



AT 



CS^^(e) = 2S(-e, 1 - e) [By,{-2E, d'^) - ^ a„ 5y„(n - 2e, d'^ 

n=l 
N 

- 2B{-E, 2-e) [By,{l - 2e, d'o) -Y,^n By, {n + 1 - 2e, < 

n=l 

In these eqs. By,(a,P) is the incomplete beta function 

ryo 
Jo 

while O*^*''^^ (cos x) in eq. (5.54) above denotes the angular integral 

n^'^''\cosx) = j d(cost?) (sini?)-2^y" d(cos^) (sint^)-^-^^ 



(5.55) 



(5.56) 



(5.57) 



X (1 — cost9) *(1 — cos X cos-!? — sin X sin COS 93) 
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This was computed in ref. [27] and we have (note the different normahzation of the angular 
integral here as compared to [27]) 



0(^''=)(cosx) = 2^-'- 



' vr 



r(l - 2e) 



B{l — i — e,l — k — e)2Fi [i,k,l — e 



1 + cos X 



(5.58) 



r2(l-e) 

Furthermore, from eq. ( 5.42| ) it is easy to see that 

cosXi^ = l-2%Q. (5.59) 
Next, consider S^f (y^^ Q;e) and CS^^(e). Using eqs. (|1|, |5^) and (|]52|) we find 



Sif (^ai.Q; i^) — '^Yai 



•'^27rr(l - 2e) 



cosx« 



(-2e, do + i) - ^ a„ (n - 2e, (io + | 

n=l 

B{-£,1 - e) - B{1 - e,l - e) 
|i?,„(-2e, d'o) + i?,o(-2e, d'^ + 1) - 2Bj,„(-2e, d^, + 



2> 



AT . 

^ a„ [Syo (n - 2e, d'o) + (n - 2e, do + 1) - 25j^o (n - 2e, d'o + i) I 

n=l ^ 



TV 



with cosx^^ = 1 - 2y„^ Q, and 

CS^^(e) = 2B{-e, 1 - e) [By,{-2e, d'^) - By,{n - 2e, d'^) 

n=l 

Finally, using eqs. (Q) and ( p3| ), for S^^(e) we obtain 

N 

Siiie) = -Bi-e,-e)\By,i-2e,d'o)-^a^By,{n-2e,d'o) 



(5.60) 



(5.61) 



(5.62) 



n=l 



The choice of a^i. Before presenting the e-expansion of the soft- type integrals, let us 
consider how the a„s appearing in eq. ( |5.31| ) are fixed. The crucial observation, as discussed 
in the Introduction, is that we must obtain an integrated approximate cross section which 
has a universal collinear limit. However, the colour-connected pieces of the integrated 
approximate cross section can ruin universal collinear factorization as explained below 
eq. ( |1.6D . Now as eq. ( ^.34 ) shows, the integrated soft terms S^"^* will enter the integrated 
approximate cross section multiplying colour-connected matrix elements (see eq. (|6.17|) ). 
Thus, in order to guarantee that the integrated approximate cross section does have a 
universal collinear limit, we need to insist that 



JpJl 



(5.63) 
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for any j, I and p. This condition implies certain constraints on the a„ that appear in the 
function P{yrQ) in eq. ( ^.31 ). 

To understand these constraints, consider first the case when j, I and p are all final-state 
partons. Then to satisfy eq. ( ^.631 ) we need the following equality to hold: 



(5.64) 



As is clear from eq. ( ^.36 ), the variable Yi^^q is homogeneous in the momenta and p^, 
thus in the collinear limit, when pj zpj^ and — > (1 — z)pji, we simply find 



C~-Y-- 



^jl^pl,Q 



Y~ 
pjl,Q ' 



(5.65) 



and eq. ( 5.64| ) is trivially satisfied. The same argument shows that if now p is an initial-state 
parton, 

also holds. 

Next, let p and e.g. j be final-state partons, while I is an initial-state one. By eq. ( 5.63| ) 
we require 

C'.i^pfi'YplQ-^^) = C,7S^/(ypz,Q; e) = C.,Sf/(yzp,g; e) . (5.67) 

Of course eq. ( |5.65| ) is again satisfied (recall that if / is an initial-state parton then pi = pi), 
but it is obvious from eqs. ( ^.54 ) and ( ^.60D that S^^(y;e) and Sjp{Y;e) are in general 
different functions of the variables Y and e. We have included the P{yrQ) function in the 
definition of the soft subtraction term, eq. ( ^.27 ), precisely to have the freedom to enforce 
eq. (|5.67| ) above. In particular, we if require that {Y;e) and Sjp{Y;e) are equal, then 
by virtue of eq. ( pSD , eq. ( p7| ) holds. We find that S^f{Y;e) = Sl^{Y-e) is satisfied if 
we demand that 



N 



By^ (-2e, + 1) - ^ a„ By^^ (n - 2e, (1'^ + !) = 

n=l 

N 

= Byo (-2e, d'o + ^) - ^ fln By,^ (n - 2e, + ^) , 



(5.68) 



n=l 



and 



= By,{-2e, d'o) + By,{-2e, d'^ + 1) - 2By,{-2e, d'o + i) 

N 

an [By, (n - 2e, d'^) + By^ (n - 2e, 4 + 1) - 2By, (n - 2e, d'^ + I] 



(5.69) 



n=l 



Finally, if p and I are the two initial-state partons while j is a final-state one, eq. ( 5.63| ) 
implies that we need to have 



^jiKj0^pj,Q'' ^) ~ Cj^Spf (e) . 
In this limit we have p^ — > (1 — x)pf, so we find 



(5.70) 
(5.71) 
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because p and I are the two initial-state partons (see the discussion below eq. ( 5.36| )). But 
as can be seen from eqs. ( 5.6d| ) and ( 5.62| ), SpJ(l;e) and Sp[(e) are generally different 
functions of e. Demanding their equality given eq. ( p9D impl les 



TV 

By,{-2e, d'o + ^)-Y,an By,{n - 2e, d'^ + i) 

n=l 



N 



(5.72) 



Byg{-2e, do)-\an By^ {n - 2e, d^ 



n=l 



We have used that 



r2(i 



B{ 



Thus, we find that if we choose the o„ such that eqs. ( 5.6^ ) and ( [5.69| ) as well as 
eq. (|5.72| ) are satisfied, then 

\Y;e) = SifiY;e) (5.74) 



cFF^ 



and 



sif(i;^) 



(5.75) 



so the coUinear limit of the integrated soft counterterm is universal. 

To solve for the a„, note first that eqs. ( p.68| , p. 691 ) and ( |5.72| ) are not independent. In 
fact, eqs. ( |5.68| ) and ( 5.72| ) clearly imply eq. ( ^.69| ). Therefore we only have two independent 
equations and thus the minimal number of a„s is = 2. Here we limit ourselves to this 
minimal solution and use eq. ( ^.72|) to get 



Byg{-2e,dQ + ^) - -Bj/o(-2e,do) = By^^{l - 2e,dQ + ^) - By^^{l - 2e,dQ 

+ [By, (2 - 2e, 4 + i) - By, (2 - 2e, d'^ 

and take the difference of eqs. ( 5.72| ) and ( p.68| ) to obtain 

By,{l - 2e, d'o) = By,{2 - 2e, d'o)ai + By,{5 - 2e, 4)a2 



ai 

0-2 



(5.76) 



(5.77) 



as the two independent equations determining oi and a2- In deriving eq. ( 5.77| ) we use the 
identity 



By,{a,l3)-By,{a,f3 + l) = By,{a + l,P), 



(5.78) 



which is easy to show starting from the integral representation of the incomplete beta 
function as given in eq. ( 5.56 ). It is straightforward to solve eqs. ( 5.76 ) and ( 5.77| ) for ai 
and a2, and we give the solution in appendix |A.2| . As remarked before, oi and 02 depend 
on yoi d'^ and e. It is important to point out that for = D'q + d'^e with D'q an integer 
{D'q > 2, see appendix A of ref. [23]), ai and 02 are finite as e — > 0, which is necessary for 
eqs. (5.27)-( 5.2g| ) to make sense as subtraction terms in d = 4 dimensions. 

Finally, we remark that using eq. ( |5.77| ) to simplify eq. ( ^.551 ), we find that the second 
term of the latter equation vanishes, and thus we have (compare eqs. ( |5.55| ) and ( [5.61| )) 



CS^^(e) = CS^^(e). 



(5.79) 
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Final results. Having determined the a^, we can write the integrated soft-type countert- 
erms explicitly. Here we limit ourselves to recalling the pole parts of the combinations 

S^^(y; e) + CS^^(e) , Sif (y; s) + ^ [CS^^(e) + CS'^{s)] , and SiUe) + CS^^(e) , 

(5.80) 

as it is these combinations which finally enter the insertion operator (see sect.^]^ below). 
We find that the pole parts are independent of yo and dg and read 

Sf/(y; e) + CS^^ie) = Sif{Y; e) + ^ [CS^^(e) + CS^^(6)] = ^ In y + 0{e') (5.81) 

and 

Sii{e) + CS'^{e)=0 (5.82) 

exactly. (Notice that 2B{—e, 1 — e) = B{—e, — e), so that the right hand sides of eqs. ( |5.61| ) 
and ( 5.62| ) coincide up to the overall sign.) The finite part of eq. ( |5.81| ), with cIq = 3 — 3e, 



is given in appendix IA.4 

6. Implementation of the subtraction scheme 

The unpolarized fully differential real-emission cross section with two incoming partons a 
and b multiplied by the jet function reads 

da^f,{pa,Pb)Jm+l =M V --7 r#m+l(M;Pa +Pfe)-7TTT7T 

S~^,x'^iPa-Pb) Uj{fa)Uj{fb) 



X abi{P}'^Pa^Pb)\'^Jm+l{{p};Pa,Pb) 

^{m+1} 



Here JV denotes all QCD-independent factors, the factor l/{ui{fa)uj{fb)) accounts for the 
average over the polarizations and colours of the initial partons and $(pa 

• pb) is the flux 

factor, which fulfills the following scaling property 

^{r)Pa-Pb) = rj^{Pa-Pb)- (6.2) 

Finally, X^j^+i} denotes a sum over the different subprocesses that contribute and S^^+i} 
is the Bose symmetry factor for identical particles in the final state. 

6.1 The approximate cross section 

In order to compute the NLO cross section, we write it as in eq. ( ^.11 ), with the approximate 
cross section dcr^' {pa,Pb) given by 

<^(^^^iPa,Pb)Jm=^f V -^--^ -d4>m+l{{p};Pa+Pb) , . 

{m+1} ^g3) 



'-'{m+1} 



■^l\MZ+l,abi{P}-^Pa,Pb)\^^Jm, 
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Figure 1: Scatter plots of 100 sequences of 30 points each, demonstrating the convergence of the 
approximate cross section to the real-emission cross section in the process gg — > Ag. The ratio R is 
defined in eq. (^). We have set uq — y'^ = yo = 1 and Dq = D'^ ^ 3. From left to right the plots 
show the C34 final-final coUinear limit, the Cai initial-final coUinear limit and the S3 soft limit 
respectively. The labels of the two incoming gluons are a and b, the final-state gluons are labeled 
1 through 4. 



where we define 



+ Y.^ar{{p}-^Pa,Pb)Jra{{p}t"^\Par,Pb)+ [Sr{{p}\Pa,Pb) (6.4) 
ae/ \ 

- ^Cf/ S^{{p}'^Pa,Pb) - 'YCi^S^{{p};Pa,Pb) \jmi{p}^Z^;Pa,Pb) 

As previously, F and / denote the set of final- and initial-state partons respectively. 

We note that the singly-unresolved approximate cross section defined in eq. ( |6.3| ) above 
is fully local, i.e. all azimuthal and colour correlations are properly taken into account. This 
means in particular that we can test numerically the convergence of do"^'^ Jm to the real- 
emission cross section dcr^^Jm+i in any singly-unresolved limit. As a check of the proposed 
scheme, we have examined the process gg — > ig, where non-trivial azimuthal and colour 
correlations are both present. By generating sequences of phase space points tending to a 
particular limit, we have confirmed numerically that 



dfjR J4 



(6.5) 



in all one-parton unresolved limits. This is illustrated in fig. |l|, where the scatter plots show 
100 sequences of 30 points each, starting from random phase space points and converging 
to a given limit. One sequence of points is highlighted in each case for transparency. 
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6.2 The integrated approximate cross section 

In order to compute the integral of the approximate cross section, let us write it as follows 

<io-^b^{Pa,Pb)Jm = da^b^'{pa,Pb)Jm + da^b^" {pa, Pb) Jm + dafb^"'{pa,Pb)Jrn , (6.6) 

where the three terms correspond to the final-final collinear, initial-final collinear and soft- 
type terms respectively. 

To evaluate J^,ida^^^'{pa,Pb)Jm, we use the phase space factorization property of 



eq. (5.7) and perform the integral to find 

d<^^^'iPa,Pb)Jm= I Yj ^TT^ 

1 -^"^ im+1} ■ 



-d4>„,{{p}t'^;Pa+Pb) ^ 



x^E|:£<^)^Cr(.,,.)Tf. (6.7) 

X \M^ZbiiP}'^i:^■^P-^Pb)\'Jn^m^'^■,Pa,Pb) ■ 

This expression is not yet in the form of an m-parton contribution times a factor. In order 
to rewrite eq. ( |6.7D in such a form, we still have to perform the counting of the symmetry 
factors in going form m to m + 1 partons. This counting was performed in ref. [5] (see 
their eq. (7.19)) and here we only recall the final result 



ir GF 



f~=CI f~=q f~=a 

im ^ ' ir&F 

^ ir ^ 



(6.8) 



where the upper indices on the sums indicate that we may go from m to m -|- 1 parton 
final states by (i) adding a gluon or (ii) exchanging a gluon for a quark-antiquark pair. 



Substituting eq. ( |6.8| ) into eq. (6^) we find 

/ da^b^^{pa,Pb)Jm= ^Y. ^(r!^ ^ n #m(M^'^ Pa + Pfe) -Ty^ 
Jm+l Jm t-.'i^yPa-Pb) ^\Ja)^ 

ir eF 

X \MZ]abmi^^-^Pa,Pb)\'Jrn{{p}g'^-,Pa,Pb) , 



(6.9) 



where we have defined 



Cf (x; e) = Cf (x; e) ee Cf/(x; e) , (6.10) 



Cf (x; e) ^ icf/(x; e) + nf Cf/(x; e) . (6.11) 
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For further reference, we note that the pole parts of Cf (x; e) may be written in the unified 
form 



(6.12) 



with the usual flavour constants 



(6.13) 



Next we consider the integral , ^ '^o^h^"\Pa,Pb)Jm- Using the phase space convolu- 



tion property of eq. (5.19) and performing the integral we find 



m+l 



Jm ( I 1 W 



dx 



'i<Pm{{p}';Z^;xpa+Pb) 



{m+l} ■ 



1 



^{Xpa-Pb) 



(ar) V''^) ^) ^ ar 



(ar). 



(6.14) 



To get eq. ( |6.14 ) in the form of an m-parton contribution times a factor we again have to 
perform the counting of symmetry factors. This counting was performed in ref. [5] (see 
their eq. (8.18)) and is anyway trivial in this case, 



{^+1} ^{m+l} 



fr {m} ^ ' 



(6.15) 



Using the flavour sum rules, we can rewrite the sum over the flavours of r as a sum over 
the flavours of the initial-state parton ar (c = far), and eq. ( |6.14D becomes 

^^^li^''{Pa,Pl^Jm = j^MY^ j dx ^ _ Yl ^(l>m{{p}''^;xpa + Pb) 

m+l m {m} ^ c 



1 tts 



(^2 \ ^ 
Q2 ) ^acK^'-i ^) 



(6.16) 



X \M';^\f,{{p}';ii^;xpa,Pb)\^Jm{{p}ia^;xpa,Pb) + {a ^ b) . 



Finally we turn to the evaluation of f^^^i da^^^"\pa,Pb)Jm- Inserting the phase space 
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factorization of eq. ( |5.32|) and performing the integral we obtain 



/ R A 

/ dcr^fe ""{Pa,Pb)Jv 
Jm+1 



Jm r ,1 ^[Pa-Pb] UJ[fa)U. 



1 Os 



{m+1} 



2 \ £ 



E 



+ 2Y,Si[{Y^:,^Q;e){M^^l,m';:^-,Pa,Pb)\Tan^^^^^ 
+ 2Y,Si^{Y,j^Q;e){M^^l,{{p}^:}-,Pa,Pb)^^^^ 



Jm{{p]'>i^;Pa,Pb) ■ 



(6.17) 



Using colour conservation (eq. (^^)), we may combine the integrated soft and soft-collinear 
terms and write eq. ( |5.17D in the following form 



/ R, A 

/ dcJ^fe ""{Pa,Pb)Jn 
Jm+1 



/ E 7^ — Td(t>U{p}t^;Pa+Pb)^j^ 

Jm r ,1 ^{Pa-Pb) ^UaY 



1 as 



{m+1} 



E 



E E sf/(%Q;£)(>i2!.,({p}^;^;Pa,p.)iT,niA^^;„,({p}M;p,,p,)) 



Sei — 



(0) 



2 \ e 



+ 2^Sif(y,._Q;e)(All^;,,({p}M;p„,p,)|T,T,|M!^!,,(M(^)^^ 



+ 2j]S^f(y,._Q;.)(A^^;„,(MM;p,,p,)|nT,|A^;^;^,({p}W 



(0) 



+ 2Si^(e)(A^(^;,,(MM;p,,p,)|T„T,|A^;^;,,({p}M;^..,p,)) 



(0) 



(6.18) 



where as claimed, the integrated soft and soft-collinear counterterms appear only in the 
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following combinations: 



Sf/(y;e) = Sf/(y;e) + CS^^(e), 

Sif (y; e) = SifiY; e) + ^ [(S^^{e) + CS'^ie)] 



(6.19) 
(6.20) 
(6.21) 



Recall from sect. 5^ that in fact S^^(e) = and Sfjf'{Y;e) = Slf{Y;e). Furthermore, the 
pole parts of the latter functions are very simple: 



Sf/(y;e) = S^f(y;e) = -lny + 0(£'^). 



(6.22) 



To write eq. ( |6.18| ) in the form of an m-parton contribution times a factor, we use the 
relations between the m- and (m + l)-parton symmetry factors, Sj^+i} = {rUg + l)5|m}, 
and note that for each r of the mg + 1 final-state gluons the expression in the square 
brackets in eq. ( |5.18 ) contributes the same. Therefore, 



{m+1} ' rGF {m} ' 



and we obtain 



m+l 



Jm J. ^(Pa • Pb) ^{faY 



1 Os 



Se( — 



2 \ e 



+ 2Y,Kf{yal,Q-^e){Mf^^^,{m}-,Va,Pb)\TaT.^^^^ 
i£F 



ieF 



+ 2Si^(e)(A^j::;,,({p}^^p„p,)|T,T,|A^^;,,({p}^^p„p,)) 



lir). 



(6.23) 



<An({p}m ;Pa,m) • 
(6.24) 



Collecting eqs. dJ, |6l6t |6l2l ), adding eq. (U), and using colour conservation, we find 
that the sum of the integrated approximate cross section and the collinear counterterm can 
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be written in the form 



Jm+l J m 



Jm f 1 ^\m} Jo 



1 rl 



dx / dy 



^{xypa-pb) 



{{P}m'^xpa + ypb) 



1 



^{fc)uj{fd) 



X ^^nlcd^{P}m^ ^Pa, yPb)\I''^'^'^{{p}rn;Pa,X;Pb, W, S; f^F)\^T,cdi{P}m, XPa, VPb)) ■ 

c,d&I 

(6.25) 

The insertion operator 1°-^''^'^ depends on the colour charges, momenta and flavours of the 
QCD partons. Its explicit expression can be written as follows 

(2 \ ^ 



+ ( - Cf {y^Q■,e) + 2Sif {Ya^,Q■,e)jTaT, 
+ {-Cfiy,Q;e) + 2SifiYh.,Q;e))TkT, 

+ 2SiiiE)TaTb 



+ 



SachdS{l - x)5{l - y) 
+ [-CUx;e)-j;^(^^^P-%x)-K-U^. 



(6.26) 



TcTt 



'^6d<5(l - y) 



+ 



. i&F 
+ 



TdTi 



Ti\e 



P'%y)-Kli{y) 



TdTa 



5ac5{l - X) 



It is not difficult to check that the pole part of 1°-^''^'^ correctly cancels the pole part of 
the virtual cross section. Using S^^(e:) = and substituting the following identities into 
eq. iK2 



(- cf {y^Q;e) + Sf/{Y,k,Q;e))T,Tk = ^ 



i,keF 
i^k 



1 fTf 1 , 



TiTk + 0(e°) , 



(6.27) 



-Cf(y,Q;e)+2Sif(y„,,Q;e 



1 



(§- + 7^) y^-^ ^ (2/a/-l)+0(e°) , (6.28) 
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7^{-i + lcz] Sac 5(1 - x) + 0(eO) , (6.29) 



we find that the pole structure of eq. ( |6.26| ) is exactly as in eqs. ( |1.3D and ( |1.4[ ). Therefore the 
sum da^j^^{pa,Pb)+da^i^{pa,Pb] IJ-p) cancels all the singularities in the virtual contribution 
da^f^{pa,Pb)- The cancellation of poles is a strong check on the correctness of the proposed 
scheme. 



7. Conclusions 

Extending any NLO subtraction algorithm to NNLO accuracy in a process-independent 
way is non-trivial because the integrated singly-unresolved approximate cross section may 
fail to have a universal collinear limit. The subtraction scheme in ref. [20] leads to an 
approximate cross section that does not suffer from this problem and indeed that scheme 
can be included as part of an NNLO subtraction algorithm without any changes [16, 17]. 
However, that scheme is defined only for processes with no coloured particles in the initial 
state. 

In this paper, we have presented the generalization of the subtraction algorithm of 
ref. [20] to processes with hadronic initial states. By matching the known factorization 
formulae for the collinear and soft limits of QCD squared matrix elements and by carefully 
extending the matched expression over the full phase space, we have explicitly defined 
a singly-unresolved approximate cross section do"^'^ that is completely general (process- 
and observable-independent) and fully local (i.e. all azimuthal and colour correlations are 
properly taken into account). Furthermore, the integrated approximate cross section obeys 
universal factorization properties in the collinear and soft limits. It is then possible to build 
a singly-unresolved approximation to the integrated approximate cross section do"^'^ in 
a process-independent fashion, which is necessary for regularizing the real-virtual contri- 
bution in an NNLO computation. 

We emphasize that in order to define an approximate cross section whose integrated 
form obeys universal factorization in the singly-unresolved IR limits, we had to consider a 
'non-minimal' extension of the (soft and soft-collinear) limit formulae over the phase space. 

We have performed the integration of the approximate cross section and found that 
its pole structure is exactly the same, but with opposite sign, as that of the one-loop 
squared matrix element. Thus the integrated approximate cross section correctly cancels 
the poles of the virtual correction. As a further check, we have examined numerically the 
local convergence of da^^^Jm to the real-radiation cross section, da^^Jm+i, in the process 
gg Ag and found that their ratio tends to unity in any singly-unresolved limit. 

All analytic formulae, relevant for constructing a numerical program to compute cross 
sections at NLO accuracy using the new subtraction algorithm have been presented, and 
we anticipate that any such implementation will be usable as part of an NNLO calculation 
without any modifications, once the scheme is defined at NNLO accuracy. Setting up an 
extension of the present algorithm at NNLO seems straightforward conceptually, but it 
nevertheless poses a major technical challenge and is therefore left for later work. 
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A. The integrated subtraction terms 



In this appendix we give the finite parts of the various integrated subtraction terms. We 
set the exponents do and dg to do = dg = 3 — 3s. We choose to consider this particular 
value for the reasons explained in app. A of ref. [23]. 



A.l The integrated collinear subtraction terms 



A. 1.1 Final-final collinear 

We denote the finite, 0(£°) part of the integrated final-final collinear subtraction terms as 
nnCf/{x;e). For do = 3 - 3£ we find 



flag 5ao 
12 



;(x-i) 

7ao 



3c4 3a3 _ 

8 2{a;-l) 3(x-l)2 



nag 



1 5ao I 2ao i 13ao i ( o i 



6(x-l) 6(i;-l)^ ' 6{x-iy-^ ' 3{x-iy 



13a§ 13a§ 
6(a;-l) 12(a;-l)2 " 

3 1 



2(a;-l) 3(i;-l)^ 



^ ^ + of.r-Ui + ^(^-U^)H{0;ao) + ( - T ~ + 3{x-l)'^ + 3(x-l)a ~ 2{x-l)i' 

4 q4 r|„,3 9„,3 



3(x-iy ^ 2(i-l)4 ^ 3(3;-l) 

3(^)ff(0;x) + - 2)H{0-ao)H{l;x) + 

v3 



a, 







I- ^ -I- 

2{x-l) ^ 2 ^ 



2a; 







2a: 







8a; 



3ao 2ag «n i «^2 , 2a(, 2ao , 2an _ 2an , 



2a^ 



rr-1 3(a;-l)^ ' 
3 



3{x-iy^ 3(a:-l)a + 2(x-l)4 3(a:-l)5 

2)i7(0,ci(ao);x) + (2 - ^-A^)if(l,0;a;) + 



+ 



+ f )i?(ci(ao);a;) +4iJ(0,0;a;) + ((^■ 
2)F(l,ci(ao);^)- 



(s-i)^ 



2.g(ci(ao),ci(ao);3:) _ tt^ 

(a:-l)5 3(a;-l)5 



^ + 



2 ' 



(A.l) 
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and 



J^inCq^ {x; e) — ^ 



^0 , 2ag lOal, _ Sag «g 

3(xao-2ao-a;) 3(x-2) 3(xao-2oo-2:) 2(x-2) 6(a;-l) 



20a;^ 



5ai 



3{xao-2ao-x) ^ 3{x-2) 



_i — _|_ lia _|_ 



0—2. 



^""0 ^Qq ^""0 _|_ ^"0 I 

3{a;-2) 3(a:-l) 3{a;oo-2ao-2^) 9(x-2)2 9(x-l)2 

-2 



20oS 



2q3 



9(x-2)3 



+ 



10a2_ 



+ 



3(x-2) ^ a;-l 3(a;ao-2ao-3:^) 3(x-2)2 3(x-l)2 3(a;-2)» ^ 3{x-l) 



40a, 



3(x-2) 



2ao 



3{x-l) + 3{xao-2ao-x) ~^ 3(x-l)2 3{x-l)» ^ 3(x-l)4 ^ 3(x-2)5 



2ao 



320 



3(x-2)'* 

320 
3(x-2)6 



J- 9^2 
Tp" + I + 3{a;-2)5 3(x-2)e 

g(o; x) + " 3-(;-7)r- + 



i_ 2ao 2ao , 2ao , 160ao _ 8oo i 



160 _|_ 320 



2H(ci(ao);x) 



2 ^2 



160 



3{a;-l)5 ^ 3 3{a;-2)5 



160 I 320 



3(x-2)'' 3{a;-2)B 



i?(c2(ao);a;) + 



160 



3(a;-2) 



17 + 



ln(2) - f 



(A.2) 



The other two integrated colhnear functions are not independent of eqs. (A.l) and ( [A.2| ). 
Firstly, we have the trivial relationship that Cg^{x;e) and Cg^{x;£) are equal. Secondly, 
the integrated gluon-gluon splitting function satisfies 



C^f{x;e) = 2Cf/(x;e) - (1 - e)^Cf/(x;e) 
and so for the finite part we have 



(A.3) 



(A.4) 



The H functions appearing in eqs. (|0|) and (^) denote one- and two-dimensional har- 
monic polylogarithms [28,29]. Those involving ci(ao) or 02(00) were defined in ref. [24] by 
an extension of the standard basis for 2dHPL 's (see appendix ^ below). However, all HPL 's 
and 2dHPL's that appear above can be written in terms of logarithms and dilogarithms. 
We present the explicit expressions in appendix |B|. 

For oo = 1, the expressions above simplify quite significantly, and we find 



J^inCql" {x; e; ao = 1) 



+ 



+ 



24(a;-l) 12(x-l)2 ' 3(2:-l)^ ' 3(a;-l) 



+ [-'4 



2ix-l) 



+ 



3(a:~l)2 + 3(x-l)3 2(a:-l)4 3(x-l)= 



and 



^-3(iV Mx)+21n^(x)+ 2-^ 



(a;-l)5 



Li2(l-x)-^ + f , 



(A.5) 



J^mCf/(x;e;ao = 1) = ^ 



1 -^ 5 , 2 



20 



1 I 40 I 2 I 160 I 



3{a;-l)3 3(a;-2)4 3(a:-l)4 3(x-2)5 \^3(x-2)5 3(a:-2)6 

ln(x) — 



6(x-2) 6(x-l) 9(x-2)^ 9(a;-l)^ 9(x-2)^ 

ln(2) + ' 2 



160 _|_ 320 



I 2_ 

3(a:-l)5 ^ 3 



160 _ 320 

3{x-2)^ 3{x-2)'^ 



(A.6) 
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Finally we point out that despite the appearance of factors of 1/(1 — x) in eqs. (AJ, A.2 



A. 5) and ( |A.6| ) above, the integrated collinear functions are finite at x = 1 as expected, 
and we have explicitly 



l;e) 



19a| _ 9a4 ^ _ 13^2 3^ 
50 4 9 2 2 



5 2 3 



In(ao) - In^(ao) 



2 ' 



and 



34a; 



^ + Hag - 36ao + 



+ 



3(ao+l) 



202 
9 



For both ckq = 1 and x = 1, we have simply 



^mCf/(x = 1; e; ao = 1) = If - ^ , 
HnCl[{x = 1; e; ao = 1) = ^ 



1786 _^ 1601n(2) 



45 



2aQ + 6ao 
(A.7) 

+ if ln(ao + l) 

(A.8) 

(A.9) 
(A.IO) 



In passing we note that the forms of these functions as given in eqs. ( [A.l| ) and ( [A .21 ) (or 
in eqs. ( [A.sp and ( [A.(^ ) for ao = 1) are not particularly well suited for direct numerical 
evaluation very close to x = 1 because of issues of numerical stability. However, since the 
functions are smooth at x = 1, it is straightforward to develop simple approximations to 
any desired accuracy around this one point, e.g. most simply by Taylor-expanding around 
X = 1, with the leading terms given in eqs. ( |A.7| ) and ( [A.8| ) (respectively in eqs. ( [A. 9]) 
and ( [Al0|) for ao = 1). 

A. 1.2 Initial-final collinear 

The finite part of the integrated initial-final collinear subtraction terms is denoted by 
.HnC^^(x; e). For generic y'^, we find 



HnClg{x;£) =— 



HnClq{x]e) =^<j [x^ + (1 



Cf \l + [l-xf 



ln(l -x)+ ln(l - x)e(x - (1 - y'o)) 



+ lny^e((l-y^)-x) 



(A.ll) 



ln(l -x)+ ln(l - x)e(x - (1 - yo)) 
+ lny^e((l-y^,) -x)l + 2x(l - x) 



(A.12) 



TinC^ Jx] e) =2 



ln(l — x) 
1 — X 

{(1 + x) 



+ 2 



Infl — xl 



+ ^e((i-y^)-x) 



ln(l - x) + ln(l - x)e(x - (1 - y'^)) + lny^e((l - y'^) - x) 
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FinC^ {x; e) =2 



- 1 + X 

/ ln(l -x) 
V l-x 

l-x 



} + <5(l-x) 



vr 
"6 



+ ln'yo) , 



(A.13) 



+ 2 



ln(l — x) 
1 — X 



-y'o 



+ ^e((l-y[,)-x) 



+ 2 



1 + x(l - x) ) ln(l - x) + ln(l - x)e(x - (1 - v'q)) 

^2 



+ lny^e((l -y'o)-x) +5{l-x)(-^ + In^ y'^ 



(A.14) 



Above the '1 — t/q' prescription is defined by its action on a generic test function g{x) as 
follows 



J dx5r(x)[/(x)]i_^/ = ^ ^ dx [g{x) - g{l)]f{x) 



(A.15) 



These integrals were first computed in ref. [25]. However, our notation is sufficiently 
different form the one employed in [25] that a direct comparison is not entirely straight- 
forward. Therefore, we give the expression for C^^(x;e) in terms of the various functions 
introduced in ref. [25]. We find 

CUx; e) = -^P"^x) + K-'ix, y'o) + K"\x, y^) + P^'{x) In x + 6-'6{l - x)Va{y'o, e) + 0{e'), 
^ (A.16) 
where P"*(x) are the four-dimensional Altarelli-Parisi probabilities given in eqs. ( |2.6D ~ (|2.9| ), 
while K"-^{x,yQ), K {x^y'^) and Vaiy'^^s) are all defined in [25] (note also that our y^ is 
the a of ref. [25]). 

Finally, we remind the reader that for y^ = 1, the C^(j(x;e) functions become identical 
(up to a colour factor) to the ^^'^(x;^) functions of ref. [5]. The precise correspondence is 
given by 

CUx-,e-y'^ = l)Tl = V-\x;e). (A.17) 
A. 2 The integrated soft-type subtraction terms 

For completeness, we give the explicit solution of eqs. ( ^.761) and ( 5.77 ) determining oi and 
02 for general yo and d'^: 

ai = [By,{l - 2e,d'o)[By,i2 - 2e,d'o) - By,{2 - 2e,d'o + i) 

- Sj^o(3 - 2e,do) By^{-2£,d'Q) - By^{-2£,d'Q + ^) | 
X {By, (2 - 2e, d'o) [By, (2 - 2e, d'^) - By, (2 - 2e, d'o + i) 

- By,{3 - 26, d'o) [By,{l - 2e, d',) - By,{l - 2e, d'^ + 1)] , (A.18) 

02 = [By,{2 - 2e,d'^)[By,{-2e,d'^) - By,{-2e,d'^ + \) 

- By,{l - 26, d'o) [By,{l - 26, d'^) - By,{l - 26, d'o + i)] } 
X [By, (2 - 26, d'^) [By, (2 - 26, d'o) - By, (2 - 26, d'^ + \ 
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- By, (3 - 2e, 4) [By,{l - 2e, 4) - (1 - 2£, ^ + i)] } \ (A.19) 



Note that the difference Byf^{—2e, dg) — i?yo(— 2e, + 5) is finite in e. For = 3 — 3£, the 
solution reads 



ai 



945 945 



60""" 45 300 30 1890 """ 63 63 V yO^QTb 

, ^yo\ , „4_ M 



5 



240 ^ 30 



60 210 



¥) ln(2) 

13^ 5^ _ 23^ 2^\1 ^.1. 
90 ^ 27 189 ^ 63 )\ ^ ^K'^- )i 



61;/(1 , 26^ 
252 189 



63 



4725 1575 
1890 270 



(A.20) 



and 



02 



72 6 60 315 105 7 V J- i/0 1 210 210 126 63 



210 210 



65^ , 6]£g 
126 63 



105 

30 
63 



¥) + (|-¥ + 2/o^)ln(vT^+l) + (-| + ¥-yo^)M2) 

60 210 



240^ 



61^ , 2%3 
252 189 



63 



+ VI - yo 



630 



1890 270 



90 27 



23yg 
189 



+ 



(A.21) 



As advertised, the solution is finite in e, and in the expansions above we have kept only the 
0(6*^) terms, since they are enough to compute the finite part in the expansion of Sf^{Y; e) 
and Sif{Y-e). 

We note that for = 1, the above expressions for ai and 02 simplify to 



ai(yo = 1) = + 1008 ln2 + 0{e^) ~ 9.09236 + 0{e^) , 

02(7/0 = 1) = 1734-25201n2 + O(ei) ~ -12.7309 + 0(ei) . 



(A.22) 



Next, wc present the finite parts of the integrated soft terms S^^(y;e) 



denoted by HnSl^ {Y\e) and J-inS j iY;e) respectively. For the specific value of dg 
3 — 3£, we find 

I3i3yf 



1313y^ 17yg 4yl 
1890 63 63 



_ 4^ 
63 



21 

270 90 27 



! y ^ \^120 

in(y)}/{ 

5 c:,,4 9' 



15 



I M _ Mi 

30 3 



139y^ I 
135 



\^[480 40 """ 225 6300 16200 

1^_L n 7-( 8%o° I 24I3/8 173yg 526yt'; 8584y,« 761;/g 1103y^ 

63 "T V i/O \^ 113400 ^ 22680 2520 2025 14175 945 1890 

' #)in(yr^+i) + (-4 + !t-^ 

30 60 ^ 210 252 ^ 189 63 ^ v UOy gSO 
-J-)} + (4yo - yg) iny - i ln2 y -2ln{yo) InY - Li; 



30 3 



240 

., 23yg 2yg 
27 189 63 



63r 



1890^ 



i2(i-y). 

(A.23) 
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Setting Ho = 1 leads to substantial simplification: 

^nSf/(y; e; yo = I) = ^nS^f (F; e; yo = 1) 



[2516 _ 252ln2)lny + iln2y + Li2(l-y)l . ^^'^^^ 



V 15 '''^ ^2 

Finally, we remind the reader that S^^(e) = 0, and this is true to all orders in e, for any 
yo and d^. Therefore we have Sfjf{Y = l;e) = S^f {Y = 1; e) = as well, to all orders in 
e. Clearly, !PinSfif{Y; e) and ^nS^f (Y; e) as given in eq. (|A.23[) are indeed zero at y = 1. 

B. Explicit expressions for some harmonic polylogarithms 

In this appendix we collect all one- and two-dimensional harmonic polylogarithms that 
appear in the integrated final- final collinear subtraction terms (see eqs. ( [A.l| ) and ( |A.2| )), 
expressed in terms of logarithms and dilogarithms. The weight one HPL 's are 

H{0;ao) =lnao, (B.l) 

H{0;x)=\nx, (B.2) 

F(l;x) =-ln(l-x), (B.3) 

while the weight two HPL 's read 

H{0,0;x) =^ln^x, (B.4) 

H{l,0;x)=U2{l-x)-^. (B.5) 

o 

The functions involving ci(ao) or C2(ao) were defined by an extension of the standard basis 
for 2dHPL 's in ref . [24] . The new basis functions read 

f{ci{ao);x) = — — -, f{c2{ao);x) = — — -, (B.6) 

x-ci(ao) x-C2{ao) 

with 

ci(ao) = "° -I , C2(ao) = — ^ . (B.7) 
The 2dHPL's involving ci(ao) are 

H{ci{ao);x) =ln(l + ^—^x] , (B.8) 

ii'(0,ci(ao);x) =-Li2 f-i^xV (B.9) 

H{l,ci{ao);x) = - ln[(l - ao)(l - x)] In ( 1 + ^ ~ | (B.IO) 

-FLi2(ao) -Li2[l-F(l-ao)a;], (B.ll) 

i/(ci(ao),ci(ao);x) =iln2 (l + l^x") . (B.12) 
Finally, the single function involving C2(ao) evaluates to 

F(c2(ao); x) =ln(l + ^^^^A ■ (B.13) 
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